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Abstract 

Numerical simulations of the chiral phase transition in the (3+l)dimensional 
0(4)-model are presented. The evolutions of the chiral field follow purely 
dissipative dynamics, starting from random chirally symmetric initial config- 
urations down to the true vacuum with spontaneously broken symmetry. The 
model stabilizes topological textures which are formed together with domains 
of disoriented chiral condensate (DCC) during the roll-down phase. The clas- 
sically evolving field acts as source for the emission of pions and a mesons. The 
exponents of power laws for the growth of angular correlations and for emis- 
sion rates are extracted. Fluctuations in the abundance ratios for neutral and 
charged pions are compared with those for uncorrelated sources as potential 
signature for the chiral phase transition after heavy- ion collisions. It is found 
that the presence of stabilizing textures (baryons and antibaryons) prevents 
sufficiently rapid growth of DCC-domain size, so observability of anomalous 
tails in the abundance ratios is unlikely. However, the transient formation 
of growing DCC domains causes sizable broadening of the distributions as 
compared to the statistical widths of generic sources. 
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I. INTRODUCTION 



There is general consent that the energy densities achieved in relativistic heavy-ion col- 
lision experiments performed or planned at RHIC or LHC should be sufficient to drive 
hadronic matter through the QCD transition into a new phase with restored chiral sym- 
metry and, probably, deconfined color. A commonly accepted scenario |L| assumes that 
immediately after (central) collisions most of the baryon number is still concentrated in 
Lorentz-contracted pancakes receding from the collision point with approximately the speed 
of light (in beam direction), leaving behind in the cm. system a rapidly expanding cylin- 
drical or spherical fireball of highly excited "vacuum". Although almost void of net baryon 
number, the energy density initially deposited in this 'fireball' is sufficient to produce during 
the subsequent cooling phase a large number of hadronic particles, mostly pions, but also 
heavier mesons and baryon-antibaryon pairs. If inside the initial hot fireball chiral symmetry 
indeed has been restored, then during the following relaxation process the relevant order pa- 
rameter, the chiral condensate, must re-increase to finally approach the nonzero value which 
characterizes the true vacuum. The time scale for this evolution of the order parameter is 
set by typical relaxation times r of the strongly interacting system (e.g. through dissipa- 
tion of energy by particle emission) which need not necessarily be closely tied to kinematic 
time scales r c of the collision and expansion time of the fireball. If it is possible to replace 
during this relaxation process the complex dynamics of the strongly interacting fermionic 
and bosonic elementary degrees of freedom by effective dynamics for an order-parameter 
field, then the time scale for changes in the relevant effective potential can be quite distinct 
from the time scales of dissipative processes which characterize the relaxation of the order 
parameter towards low-energy configurations in the changing effective potential. So, within 
that concept we can study evolutions which proceed through sequences of configurations 
which are far from equilibration. Especially, for r c <C r, (the so-called "sudden quench"), 
the system moves through maximally non-equilibrated sequences. 

The standard tool for an effective description of hadronic physics proceeding on energy 
scales of less than 1 GeV is the chiral 0(4) model. The isoscalar quark condensate (qq) (with 
finite value (qq) in the true vacuum) is combined with a corresponding isovector condensate 
(qrq) to form a four-component "chiral" field 3>, considered as space-time-dependent order- 
parameter field, and subjected to an effective action which governs the "chiral dynamics". 
(We do not discuss here the natural extension to the full SU(3) flavor group.) Fluctuations 
of $ describe the pseudoscalar pions and scalar a- mesons. In its simplest form, the "linear 
cr-model", the effective action comprises the kinetic two-gradient term and an appropriate 
time- or temperature-dependent $ 4 potential. In this form the model has been widely 
used for investigating features of the chiral phase transition, both in classical dynamics [0 , 
including dissipation and noise 0, and in quantum field theoretical approaches [[§]. 

The generic structure of this model applies to a wide variety of physical systems with 
spontaneously broken symmetry, such that particular phenomena may be expected to occur 
irrespective of the nature of the underlying physical degrees of freedom. Depending on 
the number of field components and spatial dimensions the formation of ordered domains, 
separated by topological defects or textures characterizes ordering evolutions in all such 
systems |J. This analogy has prompted the idea to use pions emitted from differently 
oriented ordered domains as indicator for the chiral phase transition itself |J . It is expected 
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that the temporary misalignment of the chiral field in finite spatial domains (as compared to 
the orientation of $ in the true vacuum) leads to anomalous fluctuations in the multiplicities 
of emitted charged or neutral pions which may provide a signature for the transient existence 
of disoriented chiral condensate (DCC). 

On the other hand, for temperature T=0 the interpretation of topologically nontrivial 
^-configurations as baryons |7|] (in linear or nonlinear realizations of the 0(4)-symmetry) 
has provided surprisingly successful models for baryon structure and dynamics. Their sta- 
bilization requires additional terms in a gradient expansion of the effective action. Their 
form is fixed by chiral symmetry, the magnitude of the few relevant low-energy constants is 
extracted from experimental information. These higher-order gradient terms incorporate at 
least partly the influence of vector mesons. Naturally, the well-established Skyrme term is 
the appropriate and simplest tool for this purpose. 

As it is well known that stable textures play a most important role in the ordering 
evolutions of multi-component fields, it appears essential to allow for their formation and 
stabilization in the dynamics of evolving field configurations. At the same time this opens 
up the possibility to describe creation and annihilation of baryonic structures during the 
evolution of the fireballs. These structures have to be supplied with a definite spatial length 
scale which is provided by the stabilizing mechanism in the effective action, which prevents 
them from collapse or indefinite growth. Using all four components of the 0(4)-field as 
independent variables implies that in the phase with restored symmetry the absolute value 
of the condensate can be very small locally, as compared to its true vacuum value. This is 
different from the nonlinear version of the 0(4)-model where $ is restricted to the 3-sphere 
$ 2 = f% and symmetry restoration can only occur through spatial averaging over different 
orientations. 

The true vacuum which surrounds the fireball sets the boundary condition for the evo- 
lution of the ordering field such that $(cc, t) approaches a fixed vector for r — > oo. In 
that case, all configurations of the nonlinear model fall into classes distinguished by integer 
winding numbers B. The same is true for configurations <&(x,t) of the linear 0(4)-model, 
as long as vanishing field vectors are excluded. While B in the nonlinear model is topo- 
logically conserved in time, in the linear version integer jumps in B may occur if at some 
point in spacetime $(a?,t) moves through $ = 0. Topologically nontrivial configurations 
with nonvanishing $ embedded in the surrounding vacuum field then are interpreted as 
baryons with baryon number B. In the nonlinear version of the 0(4)-model these 'baryons' 
(with baryon number B and definite size) are topologically stable against unwinding. In the 
linear 0(4)-model they represent metastable configurations which occasionally may unwind, 
especially if chiral symmetry is broken explicitly through a non-vanishing pion mass. As 
topological arguments do not apply to lattice implementations, B- conservation if desired 
has to be imposed in both cases (linear or nonlinear 0(4)-model) by an optional additional 
constraint, which rejects occasional configuration updates that would imply a jump in B. 

Numerical simulations based on such very specific effective actions which incorporate 
definite scales will provide answers for physical observables which characterize the specific 
physical system under consideration. So, it is not so much our aim to look for universal 
features of a certain class of models, but rather to obtain statements about specific physical 
situations which hopefully comprise the essentials of a selected class of experiments. 

If, during the initial stages of the fireball evolution, color indeed is deconfined, then 
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evidently, the applicability of the 0(4)-model is restricted to a later phase where the color- 
degrees of freedom are re-frozen into the colorless order-parameter field <fr. We assume that 
this local color- confining process does not lead to any preferred local direction of the chiral 
field. So, assuming that the relaxation of $ is driven by effective chiral dynamics, implies 
that color-confining transition and the restoration of chiral symmetry may be considered 
separately 

Although dissipative evolution of the order-parameter field is a classical concept, the 
field fluctuations around the evolving configurations acquire particle properties and must be 
quantized as emitted radiation. With the classical field acting as a source for the emission 
of field quanta it is expected that the radiation carries signatures of the source field config- 
uration. This idea underlies the search for DCC effects in the pion yields after heavy-ion 
collisions. To define the separation of the classically evolving part from the quantum fluc- 
tuations we remove all propagating terms (second time-derivatives) from the equations of 
motion for the classical configuration, i.e. the classical evolution is defined as solution of the 
purely dissipative Time-Dependent-Ginsburg-Landau (TDGL) equations. As a consequence, 
the radiation is driven by the instantaneous field velocities, from which the event-by-event 
fluctuations in abundance ratios may be obtained. 

The following discussion represents an extension of previous work dealing with corre- 
sponding features in the lower-dimensional (2+l)d 0(3) model ||. In contrast to the ex- 
ploratory nature of the previous investigations we here attempt to define the model as close 
as possible to the scales which characterize actual hadronic matter and to the physical situa- 
tion which might prevail after an ion-ion collision. In Chpt.II we formulate the equations of 
motion and the resulting emission rates. For comparison we present also the corresponding 
expressions for the abundance ratios for DCC sources and for uncorrelated stochastic sources. 
The model and the relevant parameters which will be used for the numerical evolutions are 
specified in Chpt.III. In Chpt.IV the initial and boundary conditions for these evolutions 
are defined, together with important observables measured during individual evolutions, like 
net baryon number, baryon plus antibaryon number, equal-time angular correlations as a 
measure for the size of disoriented domains. Chpt.V presents numerical results, from which 
typical exponents for growth laws are extracted. These are finally used to draw conclu- 
sions about the possibility to observe signatures of the chiral phase transition from charge 
fluctuations in the pion radiation emitted from the relaxing source. 



where the kinetic energy density T[<&] comprises all terms containing time-derivatives of 
We follow the evolution of field configurations, which at some initial time deviate from 
the global equilibrium configuration (the 'true vacuum') within some finite spatial region V. 



II. NONEQUILIBRIUM FIELD DYNAMICS AND DCC SIGNALS 



A. The TDGL equations 



We consider an effective action for an 0(4) vector field $ 




(1) 
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The equations of motion 5S/8<£> = which govern the classical evolution of the initial non- 
equilibrium configuration describe both, the approach towards the vacuum configuration in 
the interior of that spatial region and the propagation of outgoing (distorted) waves into the 
(initially undisturbed) surrounding vacuum. Although, of course, the classical equations of 
motion conserve the total energy E = J (T + 14) d 3 x, the outgoing waves carry away energy 
from the interior of the spatial region in which the total energy E = T + U initially was 
located. 

In quantum field theory the quantization of the propagating fluctuations describes mul- 
tiparticle production of pions emitted from the relaxing field configurations. In a separation 

*(aj,t) = *cj(aj,t) + <y*(aj,t) (2) 

where comprises all of the large amplitude motion of and 5$ contains only small 
fluctuations around $ c j, expansion of <S[3>] to second order in 5$ 




+ 1 f [ 5&(x',t') ( .,. , jffL „ ,/a 1 5$(x",t")d 3 x'dt'd 3 x"dt" + <D(5& i ) (3) 




2 J J K ' ' \5&(x' \t')6$(x",t" 
and variation with respect to <5$ provides the classical expression for the fluctuating part 

6*(x,t) = 5&°\x,t) + / g(x,t,x',t') ( * S ) d 3 x'dt>. (4) 

The homogeneous part 5&°\ after quantization, represents the scattering of field quanta off 
the classical configuration <& c / (i.e. these are on-shell distorted waves). The Green's function 
Q(x, t, x', t') of the operator 5 2 S /(5$>5$>)\^ ; j relates emitted (i.e on-shell) radiation to the 

source term 6S/6&\^ i. There is, however, no unique way to separate the propagating 
fluctuations 5$ from a more or less smoothly evolving classical configuration <fr c z because 
both, <fr c z and 5$, are integral parts of one and the same evolving order-parameter field 
Conclusions drawn from one part only, are subject to the arbitrariness of the chosen 
separation. In any case, in an equation of motion that separately describes the evolution 
of <fr c z we expect a dissipative term to account for the loss of energy through the emitted 
radiation: 

V**/[#„] 

Then, with eq.(^), the time derivatives $> c i of the classical field determine the amplitudes of 
the fluctuations, i.e. they serve as driving terms for the emission of field quanta. 

For sufficiently small values of the relaxation constant r the dissipative term dominates 
the time evolution of <& c ;, propagating parts get damped away and we can replace (|J) by 
the Time-Dependent Ginzburg-Landau (TDGL) equation 
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(6) 

The potential energy functional U contains no time derivatives of therefore $ c / has no 
propagating parts, i.e. it does not pick up field momentum, and it provides a slowly moving 
adiabatically evolving classical background field. 



B. Particle emission and DCC signals 

For the adiabatically evolving process we consider the time t as parameter such that we 
have 

(J = I (ST) . d " X = -- f • Sx = - E e «(*) ( 7 ) 




T 



a=l 



and may define (omitting the index 'cl') 

£«(*) = -/ *a(M) 2 ^ = — <3l(*)> (8) 

as the energy carried away per time interval by particles emitted with field velocity oriented 
in a-direction. Here "( )" denotes the lattice average over a spatial L 3 lattice. Due to the 
slow motion of the source we expect the particles to be mainly emitted with low energies, 
i.e. with low momenta k ~ and energies u a approximately given by their masses u) a 
m a . In our present context we choose the 0(4)-symmetry to be spontaneously broken 
in 4-direction (selected by the surrounding true vacuum boundary condition or by small 
explicit symmetry breaking, or by a small bias in the initial configuration). Then the 
i=l,2,3 ('isospin') components of the 0(4) chiral field constitute three 'pionic' fluctuational 
fields 5$i = iTi with small (explicitly symmetry-breaking) mass 7711,2,3 = m*, while the a- 
fluctuation 5$ 4 acquires a large mass m„ due to the spontaneously broken symmetry. One 
of the three isospin directions {% = 3, say), we identify with the neutral 7To-mesons. So we 
expect the emission rates for neutral pions n (t) and for charged pions h ch {t) 

During the early stages of the ordering evolution the spatial averages (<&?(£)) wm be sim- 
ilar for all three isospin directions, but at late times with the formation of larger disoriented 
domains they might differ appreciably. The relative abundance for neutral pions 

= ■ m < 10 > 

n (t) + n ch (t) 

is free of unknown constants and could serve as DCC signature, if it is possible to separate 
in each individual event the small number of late time "DCC" -pions from the background 
of those produced during earlier stages, from decaying c's, and from other sources. 
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If angular ordering has been achieved within certain spatial domains such that the angular 
gradients are small within these domains then the right hand side of the TDGL-equation (f|) 
may be dominated by a nonlinear term <frF[$ 2 ] which drives the growth of the condensate 
without changing its direction. During this roll-down phase, we may approximate in the 
emission rates in Eq.(||) the velocities by $j(t)F[($ 2 )]. In the abundance ratios ( [T0|) 
the function F drops out, and we find in this approximation for the contribution of one 
disordered domain to the multiplicity ratio for neutral pions 

aw = -A < u > 

i=l,2,3 % 

where $ is the classical field within that domain. This result commonly is obtained from 
the coherent state formalism, and from the foregoing we see under which conditions and 
limitations it might apply to the actual dynamical process. In the numerical simulations we 
obtain the abundance ratios from Eqs.@,(P) and flT0|). However, it is of interest to compare 
the results with the consequences of the approximation (|TTD . 

From (O) the expected signals may be derived for the ideal case, where v disoriented 
domains of equal size have been formed. For one single disoriented domain, with the chiral 
field $ uniformly aligned into some direction, the pionic field components {% = 1,2,3) are 
parametrized as 

$i=l,2,3 ( sm $ cos 0; sm # sm 0, cos &) ■ (12) 

Then, within the approximation (|TT|), we find for the fraction of neutral pions relative to 
the number of all pions emitted 

/(tt ) =cos 2 #. (13) 

In an ensemble of events where all orientations of n are equally probable the ensemble 
average (/) of /(vr ) is, of course, (/) = 1/3. The probability P(f) to find in one event of 
that ensemble the fraction / of neutral pions then is obtained from 

as 

PU) = ^ (15) 

For v such domains of size V/v within a source of volume V the probability P„(f) to 
find the fraction / in a given event then is 

P v (f) = £ ... jf 1 8 (7 - l -{fi + ... + /„)) P(A)...P(U) dh...df v (16) 

with P(f) given by (p!5|). Some of these functions, for v — 1 to v — 10, are shown in Fig.|l|. 
For increasing v they approach the Gaussian distribution 



J f sm 9d9d4> = J fP(f)df (14) 
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FIG. 1. Probability distributions P v (f) for observing the fraction / of neutral pions in an 



ensemble of configurations with v disoriented domains of size V/u, according to (16). The dashed 
line indicates the Gaussian (|17[) for u = 10. 



Pu{f) 



1 



271(T Idcc 



exp 



2lJ uDCC , 



with width 



45z/' 



(17) 



(18) 



where "()" denotes the expectation value obtained with the single-domain probability (|l"5|). 
Already for v ~ 10 domains, the distributions flH3| ) are very close to the Gaussian (|I7]), with 
the maximum shifted only slightly to smaller values of /. 

For comparison, for uncorrelated sources, the probability to find the fraction / of neutral 
pions in an ensemble of events each of which emits (the same) total number n of pions, with 
equal probability p = 1/3 for each isospin component, is given by the binomial distribution. 
Its large-n limit is a Gaussian 



Pn(f) 



2ixo: 



exp 




with width 



p(l -p) 



n 



9n 



(19) 



(20) 



(The multiplicity distributions for neutral or charged pions commonly presented in experi- 
mental analyses contain also the event-by-event fluctuations in the total number n of pions. 
For the discussion of DCC effects it is crucial to look at the event-by-event distribution of 
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ratios, like /). The width (p0|) of the non-DCC distribution decreases with 1/y/n, while 
the width ( (Eg) for z/DCC events seems to be independent of the total number n of pions 
emitted per event. However, both, the number v = V/vo of DCC domains (with fixed 
average domain volume vd) present within the volume V, and the total number n of pions 
emitted from a random source (with fixed emission density T = n/V), are proportional to 
the total volume V of the fireball. So, in the Gaussian limit, v will be proportional to n 
and we cannot distinguish between (|Tj]) and (|H5|). By comparison we find n = (5/2) z/, or 
vd = 5/(2r). This means that a large fireball consisting of many DCC domains with radii 
Rd ~ |(5/(2r)) 1 / 3 is equivalent to a random source with emission density Y. For a typical 
value T = 1 fm -3 this radius is about Rn ~ 0.7 fm. Domains growing beyond that value (for 
fixed V) lead to a broadening of P(f) according to fll8D . This broadening may occur in the 
Gaussian limit (especially for large fireballs), where the number v of DCC domains present 
is too large to observe the strong anomalies in the shape of P v (f) for very small numbers of 
v. 



III. THE MODEL 

The 3+1 dimensional 0(4) model is defined in terms of the 4-component field $ = $<E» 
with <E> • <l> = 1, and the modulus field $ of mass dimension 1, with the following lagrangian 
density in 3 + 1 dimensions 

C = C(2) + £(4) + £(0) (21) 

which comprises the kinetic term of the linear a model 



£ (2) = -d^d^, (22) 



the four-derivative Skyrme term 

£ (4 ) = ~ [(d^d^) 2 - (d^dA)(d^d u ^)] , (23) 

and the potential 

Ao) = -^(* 2 -/ 2 0n) 2 + #$4. (24) 

The transition from a chiral-symmetric "hot" phase to the "cold" phase with spontaneously 
broken chiral symmetry is driven by the coefficient f 2 which multiplies the quadratic term in 
the potential (p4]). Generally speaking, f 2 is just some input function of time (and space), 
which is negative in hot and positive in cold areas. If the transition proceeds through 
at least locally equilibrated states, then f 2 is a function of the time-dependent (local or 
global) temperature T(t). We shall formally write T for the argument of f 2 and call T the 
"temperature" , although we should keep in mind that this may be just some parametrization 
of the input function f 2 and does not really imply that thermal equilibration is being achieved 
at every point in time. 
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We include in the potential an explicitly symmetry-breaking term acting in the intrinsic 
4-direction with constant (temperature-independent) strength H . The minimum of the 
potential is located at <& m in = (0, 0, 0, / (T)), related to the coefficient f 2 (T) of the quadratic 
term through 

/2(r) = / » (r) -^f)' (25) 

Small fluctuations around this minimum orthogonal to the 4-direction carry the "7r-mass" 



m-n 



<(T) = JL. (26) 

while the fluctuations in 4-direction are characterized by the "cr-mass" m a 

ml(T)=2K*f (T)+ml(T). (27) 

Conventionally, we define the "bag field" ip by normalizing the modulus of 3> to its (T=0) 
vacuum value fn = fo(T—*0) 

$ = U<p* (28) 



such that the bagfield equals unity in the physical vacuum. 

With the physical (T=0) values of f w =93 MeV, m n =138 MeV, with the standard Skyrme 
parameter e=4.25, the only remaining free parameter in £ is the dimensionless coupling 
constant k 2 . Through (B7I) it is related to the (T=0) cr-mass. For the range 6 < rn a / f n < 10 
we find a typical range of 18 <n 2 < 50. The spatial extension of a B=l skyrmion is (e/^) -1 , 
so we will be dealing with baryons of typical (baryonic) size of about 0.5 fm. (Not to be 
confused with their electromagnetic radii which receive sizable contributions from vector 
meson clouds.) 



IV. ROLL-DOWN AND DOMAIN GROWTH AFTER SUDDEN QUENCH 



A. Sudden quench 



As an input for the classical evolutions the detailed form of the function f 2 (T) in (|24"D 
is needed. If slow cooling proceeds through a sequence of thermally equilibrated states the 
temperature dependence of the condensate fo(T) may be obtained from thermal field theory, 
from loops at the bosonic or the underlying quark level. Typically, for temperatures above 
a critical temperature T c , the condensate fo(T) slowly approaches zero f (T/T c > 1) — > 0. 
According to (^5}|27]), the 7T- and a-masses then increase and come close to each other, 
m^jm a —>■ 1 for T/T c > 1, while the second-order coefficient f 2 (T/T c > 1) turns negative 
with increasing absolute value (cf. Fig.Q). Although all of this looks very much like the 
standard picture of a second-order phase transition, it probably should not be trusted too 
much near and above T c , because the chiral phase transition may be closely related to the 
color-deconfining transition, so number and nature of the relevant degrees of freedom may 
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FIG. 2. Typical features Q of the condensate fo(T), the coefficient of the quadratic term 
f 2 (T) in the $ 4 -potential, the pion mass m^, and the a-mass m a , as functions of the temperature 
T for the chiral phase transition. (Energies and temperature in units of /„-.) 

change drastically near T c , and our basic concept to describe the dynamics in terms of a 
hadronic order-parameter field may break down. 

However, if the cooling proceeds sufficiently fast as compared to the typical relaxation 
time t of the system we can impose a "sudden quench" where at time t=0 the system is 
prepared in some hot T > T c initial configuration, and then the "temperature" is quenched 
instantly down to T = 0. Then, for times t > the system evolves from its high-T ini- 
tial configuration, subject to the T=0 form of the potential fl24|) . Clearly, in this case no 
assumption about thermal equilibration, even more, no further information on the specific 
form of fo(T) is needed. 

B. Initial and boundary conditions 

Evolution of <fr c z proceeds through the TDGL equations of motion as they result from 
the potential part of fl21~D , (i.e. also from the Skyrme term only the spatial-derivative terms 
are kept). As far as the dissipative term originates from elimination of fluctuational modes 
which simulate a heat bath for the evolving field configuration, addition of fluctuational 
noise £ may be required. 

With the equations of motion (^) being purely first order in time derivatives it is suf- 
ficient to specify for initial conditions at t = only the field configurations themselves. 
As time-dependent dynamic fluctuations are not part of $> c i the high-temperature initial 
configuration is = in the interior of the spatial region in which the hot chiral field is 
located (apart from a small optional bias in 4-direction due to explicit symmetry breaking), 
and it takes on the true vacuum values on its boundary. The first step in the time evolution 
is therefore governed by the stochastic force £ alone. But this first step will immediately 
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create a configuration <fr c ; which is of stochastic nature itself. Therefore, choosing each of 
the four cartesian field components at each interior point of the lattice from a random Gaus- 
sian deviate, and fixing $ = (0, 0, 0, f n ) on the boundary of the lattice, provides convenient 
initial configurations with spatial correlation length less than one lattice unit. The width 
A of the Gaussian deviate reflects the initial temperature T > T c . It is not really necessary 
to specify a definite value for this temperature because after about one relaxation-time unit 
(after the sudden quench) the details of the initial configuration are lost anyway and evolu- 
tions proceed quite similarly, for widths chosen in the range 0.1 < A/f n < 0.5. Numerically, 
very small values of $ (i.e. small values of the width A) require very small timesteps at the 
beginning, therefore for convenience we generally choose A/f 7T = 0.3. The (T = 0) vacuum 
boundary conditions 

*(*) = (0,0,0,/*) (29) 

are kept fixed for all times at the lattice surface. This allows compactification of 3-space to 
a 3-sphere S s , and guarantees integer winding numbers. 

Starting from these initial conditions events are generated by evolving the configurations 
according to the TDGL eq.(^]). Because the time scale enters only through the <fr term, 
it is convenient to measure the time in relaxation-time units r. A large number of events 
generated in this way constitute the statistical ensemble from which ensemble averages then 
can be obtained at some point in time during the evolution. 



C. Conservation of baryon number 

The topological current is 

U = j^^ P *eabcd® a d v ® b d p $ c d°$ d (30) 

which satisfies d^j^ = 0. It allows to assign a value p(i,j, k) of the winding density p = jo 
to each point k) of the cubic lattice (or, more precisely, to each elementary lattice cube 
with lowerleft corner k)) such that the total winding number 

B= J2 p(i,j,k) (31) 

i ,j,k=Q 

summed over the whole lattice is integer and configurations can be selected with some desired 
integer value of B. We also define 

N= £ \p(i,j,k)\ (32) 

i ,j,k=0 

by summing up the absolute values of the local winding densities. Of course, N generally 
is not an integer. However, if a configuration describes a distribution of localized textures 
(and antitextures) which are sufficiently well separated from each other, then N is close to 
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an integer and counts the number of these textures (plus antitextures). In that case we can 
define the numbers N + ,N_ of "baryons" and "antibaryons" through 

B = N + -N_ N = N+ + AL . (33) 

Even if it is not integer we will in the following sometimes briefly call N the "number of 
baryons plus antibaryons" . 

In a lattice implementation the local update of k) at some timestep at some lattice 

point k) will occasionally lead to a discrete jump in the total winding number B. For 
well-developed localized structures this corresponds to unwinding textures or antitextures 
independently, such that B decreases or increases by one or more units. This will eventually 
happen even for implementations of the nonlinear 0(4)-model where the length <£> is con- 
strained to f n and B is topologically conserved, because the topological arguments based 
on continuity do not apply to the discrete lattice configurations. It can be prevented by an 
(optional) B filter which rejects such a 5-violating local field update at that lattice point 
and timestep. This eliminates all independent unwinding processes. Only simultaneous an- 
nihilation of texture and antitexture in the same time step remains possible, and, as the 
update proceeds locally at each lattice vertex it can happen only if texture and antitexture 
overlap. This 5-conserving evolution is characteristic for the nontrivial topology of the non- 
linear 0(4)-model and in this way can be implemented as an optional constraint also into 
the linear 0(4)-model. In the continuum limit it implies that a vanishing modulus $ = of 
the field vector is excluded. 

Practically, during the very early stages of an evolution, where the local gradients of the 
angular fields still are of the order of ir/a, (where a denotes the lattice constant), annihilation 
processes by far exceed local unwinding, and cause a very rapid decrease of the initially 
large number N. During that stage the definition of the global baryon number B requires 
a prescription of how to define the baryon number located on one elementary lattice cube 
(which necessarily involves some arbitrariness like mapping on geodesies). If the lattice 
constant a is chosen sufficiently small (as compared to the typical spatial size R = (ef^ 1 = 
0.5 fm of stable B = 1 baryons) such that the emerging localized extended structures 
are described with reasonable accuracy on the lattice, the small global baryon number B 
stabilizes very quickly and local unwinding no longer occurs. It turns out that R/a = 3 is 
sufficient to avoid the need for imposing an explicit constraint on B. This implies a lattice 
constant of a = 1/6 fm, which appears as a reasonable value for the lattice implementation 
of an effective low-energy model. 

As there is no universal simple scaling law (due to the ocurrence of different powers 
of gradients in £, or in other words, due to explicit scales introduced through soliton size 
and symmetry breaking) changes of the lattice constant can only partly be absorbed into 
corresponding changes in the time scale, and we have to check to which extent physical 
statements, like the physical size of oriented domains at the time when the roll-down is 
completed, are independent of the choice of the lattice constant a (for a fixed physical size 
(aL) 3 of the finite total volume). 
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D. Angular correlations 



A convenient measure for the average size Rd of (dis)oriented domains is the half- 
maximum distance R1/2 of the equal-time lattice averaged angular correlation function 



C(R) 



E E *(i,j,k)-*(l,m,n)/ E E 1 

i,j,k=0 l,m,n=0 i,j,k=0 l,m,n=0 



8<R<8+1 



8<R<8+1 



(34) 



with 5 = \J{i — I) 2 + (j - m) 2 + (k — n) 2 

Note that the definition fl34| ) contains only the unit vectors and not the full field vec- 
tors i.e. C(R) measures only the angular correlations in a given configuration. The 
half-maximum distance R1/2 is conveniently obtained by interpolating between the two 
neighbouring integer values of R where C(R) passes through C(0)/2. This lattice aver- 
aged correlation length can be obtained as function of time for any individual configuration 
(in a statistical ensemble of many events) and compared to the domain pattern of that same 
individual configuration. Analysing the domain sizes shows that Rd ~ R1/2 provides a rea- 
sonable measure for the average size Rd of (dis)oriented domains. To extract an accurate 
growth law for Rd, of course, would require to perform an additional ensemble average at 
each point in time, but for a sufficiently large lattice the spread in Ri/2{t) for different events 
is small (as long as R\n < L/2), and it is sufficient to study C(R) on an event-by-event 
basis. 

Explicit symmetry breaking (through nonvanishing H in (^4|), and through the vac- 
uum boundary conditions fl29|) imposed on <1>) causes a nonzero average ($) growing (in 
4-direction) with time, so on a finite lattice at late times it becomes numerically inconve- 
nient to extract R1/2 from C(R). Therefore, in analogy to (0), we define a correlation 
function C n (R) 

C,(R)= E E K(i,j,k)-K(l,m,n)/ J2 E 1 ( 35 ) 

i,j,k=0 l,m,n=0 i,j,k=0 l,m,n=0 

S v ' v ' 

8<R<8+1 8<R<8+1 

where 4?^ is the "pionic" (orthogonal to the 4-direction) part of 4?, renormalized to form a 3- 
dimensional unit vector 4^ • = 1. The function C n (R) measures the angular correlations 
in the pionic 3-space orthogonal to the 4-direction, unaffected by explicit symmetry breaking. 
It satisfies C n (0) = 1, and the corresponding half-maximum radius will be denoted as R v . 

The angular ordering is driven by the terms £( 2 ) and £( 4 ), which contain two, respectively 
four, angular gradients. From simple scaling arguments we expect each of these terms, 
separately, to lead to t 1//2 , respectively t 1 / 4 , power-laws for the growth of the spatial scale 
of disoriented domains. However, their relative strength introduces a definite scale into the 
system (the size (e/^) -1 of the baryons), and the presence of extended textures will cause 
deviations from the scaling growth laws. 
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V. RESULTS 



A. Roll- down times 

As a global measure for the restoration of the chiral condensate we consider the lattice 
average of the modulus of the order parameter field (if). In the interior of stable baryons 
the <p(x, t)-field will stay small (constituting the bags), while in the exterior region it will 
evolve towards the vacuum value ip=l. Depending on the size R 3 and number N of baryons 
(+ antibaryons) which finally stabilize inside the total volume (aL) 3 , the average (<p) will 
approach a value near 1 — N(R/aL) 3 . We define the typical roll-down time tfuj as that point 
in time after the sudden quench when (<p) has reached 80% of its final value. Locally, the roll- 
down towards the vacuum value (p=l in the exterior regions is driven by the potential part 
£( ) ( plQ of C. Therefore, for a large lattice, as long as the baryon(+antibaryon) density is 
small, we expect t^n to scale with k~ 2 , i.e. with the inverse strength of the <3> 4 -potential. The 
typical range for the (T=0) cx-mass ( 6 < m a / f m < 10 ) provides limits for k 2 (18 < k 2 < 50). 
This constitutes a relatively weak potential to compete with the gradient terms in C, so we 
expect deviations from a pure K _2 -power law. Similarly, the actual roll-down times are 
significantly reduced through the symmetry-breaking imposed by the finite pion mass and 
the vacuum boundary conditions which become increasingly important for smaller lattices. 
We find typical roll-down times of the order of tnD /t ~ 2-3 (for k 2 = 50 and lattice sizes 
L = 30 — 60) and £rd/t ~ 5-8 for (k 2 = 20). This shows that the potential part indeed has 
a dominating influence for the evolution of the chiral condensate. The slight prolongation 
of tnD with increasing lattice size L reflects the influence of the vacuum boundary condition 
on the growth of $ in the interior of the volume. 

B. Features of the evolutions 

Evolutions after a sudden quench show three typical distinguishable phases: 
I.) During the first phase, which extends to about t/r ~ 0.1 — 0.2 after the quench, the 
initially random configuration quickly aligns locally, such that the angular correlation lengths 
Rd ~ R n grow to about 3-4 lattice units (corresponding to about 0.5-0.7 fm). Fig. (^J) shows 
a 2-dimensional cut through a 60 3 cubic lattice with the isospin (unit)-field vectors projected 
into a plane in isospace, at times t/r = 0.02 (where local aligning is already clearly visible), 
and at t/r = 0.1, which approximately marks the end of this first phase. The condensate 
(<p) does not grow during this period, in fact, the lattice average decreases to a minimum 
slightly below its small initial value. The total baryon-plus-antibaryon number N (which 
is becoming really well-defined only towards the end of these very early stages) drops to 
values of about 100-200 on a L=60 lattice, corresponding to a (baryon+antibaryon) density 
of 0.1-0.2 fm~ 3 , while the integer baryon number B is close to or equal to zero. Already 
by the end of this short first phase the winding density has developed pronounced maxima 
and minima at the locations of the emerging (anti-)baryons, while the bag field ip drops 
smoothly from its vacuum boundary condition to its small value throughout the interior 
volume, and develops dips near the extrema of the winding density. All of these features 
are basically independent of the $ 4 -potential strength k 2 , the total volume L 3 , and the 
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FIG. 3. Typical field configurations during the early phase (I), at times t/r=0.02 and i/r=0.1 
after the sudden quench. The arrows show the isospin (unit-)vectors at the lattice points of a 
2-dimensional cut through a 60 3 cubic lattice, projected on a plane in isospace. 



symmetry breaking. The latter, however, causes a steady increase in the lattice average 
($ 4 ) of the 4-component of the vector field <fr, which initially is very close to zero. The 
energy density e = U / (aL) 3 in the random configuration at time t=0 is extremely high (of 
the order of 10-50 GeV/fm 3 ), and located almost completely in the angular gradients of the 
Skyrme term . (Note that the energy density ey = (k 2 /4)/ 4 located in the $ 4 -potential 
for $ = is only about ey=0.05 GeV/fm 3 for k 2 = 20.) By the end of this first phase the 
energy density has dropped to values of about 0.3-0.5 GeV/fm 3 , i.e. to values where we 
might hope that this effective hadronic model may become reasonable to apply. So we may 
conclude that this first phase of the evolutions rather serves to create an ensemble of initial 
conditions for the onset of the physically reasonable time-development of the model which 
begins near t/r ~ 0.1 after the sudden quench with correlation lengths of the order of 0.5 
fm. 

II.) The second ('roll-down') phase is characterized by the growth of the local ip-field 
towards unity in the spatial regions between the emerging baryons, so the lattice averaged 
condensate ((f) increases towards its final value. As discussed above, the duration of this 
roll-down phase is dominated by the strength of the $ 4 potential and by symmetry breaking. 
Fig.([|) shows a sample of ten events evolving on an L = 30 lattice with lattice constant 
a = (l/3)(e/ 7r )~ 1 = 0.167 fm, n 2 = 20, and m n = 138 MeV, with vacuum boundary 
conditions after a sudden quench. By the time of t/r = 0.1 the total 'particle' numbers 
have dropped to iV ~ 15-20 and keep falling rapidly until after some final annihilations they 
approach the fixed (net baryon) winding numbers B, which scatter between and 4 in this 
sample. The growth of the average condensate (if) is very similar for all events and differs 
only through slightly varying final values due to the different baryon numbers. 

The angular ordering increases by moving the boundaries of oriented domains, by further 
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FIG. 4. Particle numbers N, winding numbers B, correlation lengths R n (in lattice units), 
and average condensates (ip), for a sample of ten events during the roll-down phase (II) after a 
sudden quench. For comparison, the straight line (with crosses) indicates the power law 6i 0,25 ). 
The evolutions proceed on a L = 30 cubic lattice, with potential strength n 2 = 20. 

annihilations, or fusion of individual baryons into more complex multi-I? configurations. As 
long as such extended localized textures are present, they prevent straight alignment of 
the isovector part of the field in the spatial domains which surround and separate them, 
while the 4-components follow the driving symmetry breaker towards the aligned vacuum. 
Therefore, the growth of the angular correlation length is quite slow. During the roll- 
down phase it follows approximately a power growth law t a with a varying between 0.25 and 
0.3, as we might expect from the separate or at least dominant acting of the Skyrme term. 
Due to this slow growth, the average radius of aligned disoriented domains in isospace 
reaches only about 7 to 10 lattice units at roll-down time tno, and differences in tjuj due to 
different choices of k 2 do not lead to sizeable differences in R n by the end of this roll-down 
period. Depending on the baryon number remaining in the configurations after roll-down, 
the angular pionic correlation lengths R n either saturate near 10 lattice units, or (for B = 0) 
keep rising towards lattice size. However, by that time, the topologically trivial domains 
of the whole field are already oriented in 4-direction, the condensate (<&) is saturated by 
($ 4 ), i.e., the pionic components $j for 2=1,2,3 in the aligned domains practically vanish. 
In other words, outside of stable textures the field is aligned in vacuum direction before the 
(T=0) condensate is fully restored. The same holds also for the stronger potential k 2 = 50 
(i.e. m a jf lx ~ 10). 

III.) The final part of the evolution after completion of the roll-down is characterized 
by the stable textures approaching their minimal-energy configurations. Occasionally, some 
final annihilations may occur, and due to mutual attraction the locally separated baryons 
drift towards each other and finally combine into the multi-5 'nuclei' which constitute the 
classical minimal-energy configurations of the model lagrangian (|2l|) . On the relaxation- 
time scale this is an extremely slow process which extends over ~ 10 3 relaxation time units, 
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FIG. 5. Snapshots of surfaces of equal (absolute) values of winding densities during and 
beyond the roll-down phase (t/r=0.4, 0.8, 4.0, 8.0, 20, 40 from upper left to lower right) for a 
single event evolving after a sudden quench on a L=60 cubic lattice with k 2 =20. In an arbitrarily 
selected 2-dimensional cut through the lattice the isospin (unit-)vectors are also plotted to show 
the growth of aligned domains. 
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with correspondingly small field velocities. So, from the point of pion emission abundances, 
it is uninteresting. Furthermore, with propagating terms included in the classical equations 
of motion, the individual baryons or antibaryons would be able to leave the volume before 
fusing into multi-5 configurations. 

In order to convey a more pictorial view of the field evolutions we plot in figs. (|5]) 3-d views 
of the winding densities po from eq . (j30|) at six different times t/r—OA, 0.8, 4, 8, 20, 40, which 
cover the roll-down period for an individual evolution on a 60 3 lattice. The surfaces plotted 
show the equal- (positive or negative)-density surfaces, so baryons and antibaryons can be 
distinguished. Each snapshot contains also an arbitrarily selected plane in which the isospin 
field (unit-) vectors are plotted as short lines, so the growth of aligned domains is clearly 
visible. Evidently, already at the beginning of the roll-down, the individual (anti-)baryons 
are separately developed, and while the average condensate is growing towards its vacuum 
value, the angular ordering proceeds mainly through annihilation and fusion, where in this 
specific example several B = 2 torus configurations and one B = 4 cluster are created. 
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FIG. 6. Particle number N, total energy E (in units of 67r 2 /„-/e), correlation length R w (in 
lattice units), average condensate (if), average 4-component ($4/ f n ), and the abundance ratio fo 
for neutral pion emission, obtained from the field velocities according to eq. (p|-[Tc|) , for one single 
event during the roll-down phase (II). The evolution proceeds on an L = 30 cubic lattice after a 
sudden quench, with potential strength k? = 20. 



C. Pion abundance ratios 

Within our approximation scheme, the abundance of neutral pions emitted relative to 
all pions is given by the ratio of the averaged square of the field velocities (l^fLCil). In Fig.® 
we look at one typical single event evolving on a L = 30 lattice like the sample shown in 
Fig.(|j). It is characterized by baryon number B = 0. The baryon-plus-antibaryon number 
approaches the integer value N=2 near t/r ~ 1, the final annihilation process takes place 
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between t/r=5 and t/r=6. By that time the average condensate (<p) has reached 80% of its 
final value, while the average ($4//^) has almost caught up with (</?), i.e. even before the roll- 
down is finished, the field surrounding the last two remaining baryons is essentially oriented 
in 4-direction. During the roll-down phase from t/r ~ 0.2 until t/r « 6 the correlation 
length grows according to 6t ' 25 and reaches about 10 lattice units at roll-down time. 
Only after the final annihilation the growth rate accelerates. 



f ( t/T ) = f„ [T ■ d*/d t ] 




0.35 - 



fl,2,3<t/T) =fl,2, 3 [*l 

0.6 



1/3 + 0.06 ■ ( t/T ) 



1/3 - 0.06 • ( t/T ) 




FIG. 7. a) The abundance ratios /o[3>] for neutral pion emission, obtained from the field 
velocities according to eqs.(P-|To|) for a sample of 10 events evolving on an L = 30 lattice with 
k? = 20. The lower part b) shows for the same sample the ratios fi[<&] for all three isovector 
components i = 1,2,3, obtained from the field configurations according to eq.(p^). For comparison 
the growth law 1/3 ± 0.06(t/r) 3a / 2 for a = 0.25 is included. 

Figure ([]) includes the abundance ratio fo[&] for neutral pion emission, obtained from 
the field velocities as given in eqs.(||) and (|i~0l). The velocity-dependent /o[$] shows a few 
pronounced extrema, which are clearly correlated with major reordering in the topological 
structures, i.e. in this case B—B annihilations. Otherwise /o shows small fluctuations around 
the average value of 1/3. We expect that for anomalies in pionic abundance ratios the sizes 
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of aligned domains with different orientations of their isospin (l,2,3)-field components are 
important. Apart from the specific fluctuations in the velocities which occur in connection 
with annihilation or fusion processes of the emerging baryons, the pionic abundance ratios 
calculated from averaging the square of the field velocity components (0), or by averaging 
the square of the field components 

= -ifw (36) 

5=1,2,3 

show comparable deviations from the average value of 1/3. According to (|Tj^) the widths 
°vDcc °f these latter deviations (|36| ) grow with \j \fv. The number v of aligned disoriented 
domains decreases as R~ 3 oc (t")~ 3 , i.e. (J v dcc grows like ~ t 3a ^ 2 . (In this consideration 
we neglect that the effective volume available for disoriented domains is only the difference 
between the total volume and that part which is already fully aligned in 4-direction) . In 
Figs.(|7|a,b) we compare the ratios fo[&] and for a sample of events. Evidently, apart 
from the sporadic peaks in /o[$] due to B — B annihilations both distributions are reasonably 
well bound by the t 3a ^ 2 law which links them to the growing size of disoriented domains. 




FIG. 8. The rates of energy loss e a (||) (for a=l,2,3 and a = 4) through emission of pions and 
(T-mesons for a sample of 10 events evolving on an L = 30 cubic lattice after a sudden quench, with 
potential strength k 2 =20. For comparison the power-law 2.35 (t/r)~ 2-5 is indicated by the straight 
line. 

Despite this significant growth of the width of abundance ratio distributions for neutral 
pions, observable consequences are strongly suppressed by a rapid decrease of the absolute 
pion yield rates during the roll-down phase. As obtained according to eq.(H) from the 
square of the instantaneous field velocities the emission rates e a for a=l,2,3,4 are shown in 
Fig.(^) for the same sample as in Fig. (0) • These rates decrease by three to four orders of 
magnitude during the roll-down phase, whereby the explicit symmetry breaking enhances 
the time-derivatives of the 4-component in such a way that the time-integrated rates / e^dt 
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are comparable to the sum of the three pionic components Yli=i 2,3 I e idt (integrated over the 
roll-down phase). During this time interval the pionic emission rates (unaffected by explicit 
symmetry breaking) decrease approximately like ~ t _/3 with (3 ~ 2.3 — 2.5. 

The observed total number n of all pions emitted in a single event is given by the time 
integral 



n 



rco poo 

/ h(t)dt= / £i(t)/m n dt. (37) 

J to J to 2 3 



where the initial time to marks the onset of the emission. If we approximate the instantaneous 
emission rate h(t) for pions emitted from a source of volume V by a power law h(t) = 
•yV (t/to)~P, with constant emission rate density 7, we have 

n = jl—t, (38) 

If we choose for t the onset of the roll-down phase, we find for the total energy E at that 
time typically E(t /r ~ 0.1) ~ 0.3 GeV/fm 3 (La) 3 . With lattice constant a = 1/6 fm and 
lattice size L = 60 the total energy of about 300 GeV residing at time t in this (10 fm) 3 
fireball would be dissipated through emission of about 150 GeV/m n ~ 1100 low-energy pions 
and 150 GeV/m a cx-mesons. With (3 = 5/2, we obtain for the initial pion emission density 
at time to/r = 0.1, yr w 16 fm~ 3 . Evidently, most of the observed pions are emitted at the 
very beginning of the roll-down phase. 

For a random (non-DCC) source the width of the probability distribution for the fraction 
/ of neutral pions is given by the binomial result 

2 2 

a non-DCC = (39) 

where the total number n of pions emitted now is given by (|38|). For a source with u(t) DCC 
domains present at time t the width of the probability distribution for the observed fraction 
/ of all neutral pions emitted in that event is (with (|i~8D) 



1 r°° ( 4 \ 

nL h(t) {^mr (40) 

If the number v(i) of DCC-domains present at time t within the volume V is approximated 
by the power law v(t) = vq (t/to)~ 3a , where denotes the number of DCC-domains present 
at time to, we have 

aDCC = 45^7' Wlth ^ = { (3-1 ) ^ (41) 

For positive values of a the effective number of DCC-domains u e ff is reduced as compared 
to the number z/ of domains present at the time to of the onset of pion emission. (Of course, 
we have assumed here that the Gaussian approximation remains valid as long as there is 
noticeable emission strength n(t). If, however, the growth rate 3a of the average volume of 
one DCC-domain is comparable to (f3 — 1), the effective number u e ff becomes very small, 
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the Gauss approximation breaks down, and we may expect the anomalies which characterize 
P u (f) for very small values of v as shown in Fig. ([!])). 

The simulations indicate that (3 ~ 2.3 — 2.5 and a ~ 0.25 — 0.3, so we find 

Veff/vo = 0.175 - 0.5 (42) 

which implies a broadening of the distribution P(f) by a factor of 1.4 — 2.4 as compared to 
the distribution at time to- It may be noted that this result is quite sensitive to the actual 
values of a and (3: a growth rate of the angular correlation length R n ~ t 1 / 3 instead of t 1//4 , 
and the emission rates decreasing like t~ 2 instead of t~ 2 ' 5 , would lead to the limiting case 
3a = (3 — 1. Unfortunately, the result ( |4T| ) provides only a relative statement between the 
observed width obcc an d a hypothetical width at time £ , and not a comparison with the 
width ( P§D for a binomial distribution. We could enforce such a connection by assuming 
that for a = (zero growth rate for DCC domains) the source should be undistinguishable 
from a randomly emitting source. In that case we would have 



° DCC — &non-DCC* 



(3-1 



\ \(3-3a-l 



(43) 



So, the growth of the DCC domains, in competition with the decrease of the emission rates, 
still would lead to a significant broadening in the effective distribution of 7r -abundance 
ratios, as compared to the corresponding width fl3T?D for a random non-DCC source. 

The growth law for the angular correlation length R n and the roll-down times tun are 
not very sensitive to the (sufficiently large) total lattice size L, therefore the width a DC c 
of the abundance ratio distributions scales with L~ 3 / 2 . (The same holds for the binomial 
distribution). On the other hand, the amplitudes of pronounced peaks in the abundance 
ratios fo[<&] due to restructuring textures, i.e. B—B annihilation, are independent of L, while 
their number scales with L 3 . So in larger fireballs we would expect additional broadening 
of the width of fluctuations around the value of 1/3 to be increasingly due to annihilations 
of nascent baryon-antibaryon structures, without the specific anomalies which characterize 
pion emission from large disoriented domains. 



D. Slow cooling 

Up to now we have discussed the extreme case of a sudden quench where the potential 
(P3|) at time t = drops instantly to its temperature T = shape, i.e. /g(T) = f 2 for 
t > 0. Practically, however, for a smooth quench it may take a typical cooling time r c for 
/q(T) to approach its T = value f%. As long as r c is small as compared to the typical 
relaxation times r of the system, the evolutions proceed much like those for a sudden quench. 
Noticeable differences appear only if the cooling time r c becomes comparable to the typical 
roll-down time tnu following the sudden quench. 

In addition to the smooth time-dependence of fo(T(t)), the field configurations may be 
subjected to a stochastic force which describes the back-reaction of the eliminated fluctua- 
tions on the slowly evolving degrees of freedom, acting as a heat bath for the cooling fireball. 
As we have assumed in eq.(^) that the main source of energy loss is through pion and a 
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emission, with the fireball immersed in the T = vacuum, the strength of the stochastic 
force is not strictly tied to the relaxation constant r, but will be comparatively small. Let 
us therefore separately discuss the influence of the smoothly time-varying coefficient f 2 in 
the <3> 4 -potential, without any noise term. 

The first phase (I) of the ordering evolution proceeds independently from the strength of 
the <3> 4 potential, so it is also insensitive to the actual value of f 2 {T). The angular correlation 
lengths and the baryon-plus-antibaryon densities established by the end of this phase are 
independent of the cooling time r c . 

Subsequently, during the roll-down phase (II) of a slow cooling process with t c tun, 
the configurations which emerge near the end of phase (I) now have sufficient time to relax 
towards a minimum at the instantaneous temperature. As a result of the baryon radius R 
scaling as R ~ (e/o(T)) -1 , with fo(T)/ f n < 1 throughout most of the interior volume, the 
nascent textures are slightly larger, overlap and interact more strongly, and therefore anni- 
hilate or fuse more easily into multi-i? configurations. All these textures are well developed 
by the end of phase (II) although the actual condensate (<£>) is close to the momentaneous 
value of /o(T(£)), i.e. much smaller than /„.. 

During the following increase of (<&) towards f w which proceeds on the time scale r c , the 
established configurations simply follow the changing effective potential in a more or less 
equilibrated manner. The addition of a small temperature-dependent noise term causes only 
a slight retardation in the ordering evolutions. 

Altogether, the angular ordering proceeds on the time scale of the relaxation constant r, 
quite independently of the cooling time scale r c . From the point of DCC effects, the case of 
the sudden quench shows the essential features. The late-time reshaping of the established 
textures according to a changing f (T(t)) is not relevant for the observability of DCC pion 
anomalies. 

VI. CONCLUSION 

As a dynamical model for the relaxation of a hot chirally symmetric 'fireball' towards 
the cold vacuum characterized by spontaneously broken chiral symmetry with nonvanish- 
ing chiral condensate, we have investigated the classical time evolution of chiral 0(4)-field 
configurations. Starting from stochastic initial conditions, which are chirally symmetric in 
the interior of a finite spatial volume, the evolutions of the classical fields follow purely 
dissipative dynamics in their relaxation towards the physical vacuum which surrounds the 
considered volume as boundary condition. 

The model is limited by a severe restriction: The chiral 0(4)-field comprises all degrees of 
freedom of the relaxing fireball. This renders it questionable that the model can be applied 
and trusted at conditions which may prevail at temperatures close to or above the critical T c 
of the chiral phase transition. Consequently, we do not pay much attention to details of the 
initial conditions. The dissipative dynamics wipes them out, anyway, after a small fraction 
of the typical relaxation time unit. But we hope, that once the color degrees of freedom 
are (re-) confined into the hadronic chiral field, and the energy- density has dropped to less 
than lGeV/fm 3 , then the model may provide a reliable picture for the later stages of the 
evolution when the actual roll-down of the order parameter takes place. 
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While the dissipative term during the very early phase of the evolution may originate 
from the rapid expansion of the hot fireball, we assume that during the subsequent roll- 
down phase energy is lost and carried away by emission of chiral field quanta (pions and 
cr-mesons). Then the square of the instantaneous local field velocities provide a measure for 
the strength of the emitted radiation. Two features are of specific interest: The growth and 
typical size of domains with disoriented chiral field, and the possibility to observe anomalies 
in the abundance ratios of emitted pions caused by the transient formation of such domains. 

On the other hand, the model is sufficiently sophisticated to allow for creation, stabi- 
lization and mutual interactions of baryons and antibaryons in the hot hadronic gas, and 
their presence as extended topological textures is of specific importance for the ordering pro- 
cess. We consider it as one of the essentials of the present investigation to use an effective 
action which incorporates the possibility to create and stabilize these structures, because 
it is well known that textures are of crucial importance for the phenomenology of phase 
transitions. Their physical size and properties set the scale for the stabilizing and for the 
explicitly symmetry-breaking terms in the effective action. Numerical implementation on 
a 3+1 dimensional lattice then requires a lattice constant a which allows to describe these 
textures with appropriate accuracy. It turns out that a— 1/6 fm is sufficient for that purpose. 
Typical fireball diameters in heavy-ion collisions in the range from 5 to 15 fm then require 
lattice sizes of about 30 < L < 90, which are numerically reasonable to handle. 

The essentials of our results are the following: 

During the very early stages of the evolutions the angular gradient terms which drive the 
local aligning of the field vectors dominate over the $ 4 potential such that the actual roll- 
down of the modulus <3> towards its vacuum value sets in only after angular correlation lengths 
of about 0.5-0.7 fm have been established (see Fig.^c), and see the angular correlation radii 
R T in Fig.(||) near t/r=0.1-0.2). These configurations set the stage for the subsequent 
roll- down. 

During the roll-down phase the angular ordering is dominated by the higher-order gra- 
dient terms in the effective action which try to establish low-energy configurations for the 
emerging extended textures. This leads to an appreciable reduction in the growth exponent 
a of the angular correlation length R^, which measures the size of aligned disoriented do- 
mains in the isospin subspace, i.e. orthogonal to the 4-direction. Due to this slow growth, 
the correlation lengths R n reach values of only up to about 10 lattice units by the end of 
the roll-down phase (depending on the number of remaining baryons) which corresponds to 
disoriented domain radii of about 1.5 fm, independently of the total fireball size. A nonvan- 
ishing number of baryons (plus antibaryons) which by that time still may be floating within 
the developing vacuum configuration, causes the radii R^ to saturate at even smaller val- 
ues. This result, that near the end of the roll-down the disoriented domains are essentially 
pion-sized has previously been suggested by qualitatively quite different arguments [fLOfl . It 
is interesting that the purely dissipative dynamics used here, apparently leads to a similar 
conclusion. The result implies that even for the smallest fireball considered (with radius 
aL/2=2.5 fm) we expect more than at least 5 distinguishable domains to be present at the 
end of the roll-down phase. 

The realistic strength of explicit symmetry breaking in 4-direction, determined by the 
physical value of the (T=0) pion mass and by the influence of the (T=0) vacuum surrounding 
the hot finite volume, is sufficiently strong to effectuate saturation of the order parameter ($) 
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by the average 4-component ($4) even before the order parameter has reached its vacuum 
value. In other words, alignment in vacuum direction is completed before the roll-down 
is complete. This not only holds for the lower limit of the cx-mass rriu/f^ m 6 ( weak 
$ 4 -potential), but still is basically true also for m a / ~ 10 (strong $ 4 -potential). 

The number N of baryons plus antibaryons still present near the end of the roll-down 
is quite small, corresponding to densities of about 0.02/fm 3 , with large event-by-event fluc- 
tuations. Their number further decreases after completion of the roll-down due to later 
annihilations. Each one carries a typical energy of Eb = Gir 2 f n /e ~ 1.3 GeV. So, during 
the roll-down phase the energy density drops from about 0.5 GeV/fm 3 to almost zero. In 
our picture all of that is carried away by low-energy pions and cr's with the emission rates 
e a driven by the squares of the field velocities These emission rates drop rapidly during 
roll-down like ~ t~@ with (3 ~ 2.3 — 2.5 such that near t/r ~ 1, already 90% of the roll-down 
pions have been emitted. At that time, the radii of disoriented aligned domains have grown 
to at most 1 fm (6 lattice units), and the three pionic rates are still degenerate with good 
accuracy. In accordance with the growth of the angular correlation lengths like ~ t a with 
a ~ 0.25 — 0.3 the variance in the fluctuations of the abundancy ratio for neutral pions 
(obtained from the field velocities) grows like ~ t 3a . Later on during the evolution occa- 
sional strong deviations from degeneracy are due to changes in the structure of the emerging 
textures, rather than to reorientations of aligned domains. 

For the sizes of fireballs considered here (with radii larger than about 5 fm), during 
that early phase of the roll-down when most of the pions are emitted, the number of DCC 
domains present is sufficiently large to prevent significant deviations from the shape of a 
Gaussian distribution in the charge fluctuations. However, their transient formation during 
the chiral phase transition would lead to a broadening of the width of this distribution by a 



factor of y (f3 — l)/(/3 — 3cc — 1) as compared to a randomly emitting uncorrelated source. 
For the values of a and f3 extracted from our simulations this factor amounts to at least 
~ v2. Slightly smaller values for (3 and larger values for a which might still be compatible 
with the simulations could increase this factor appreciably. The slow growth of the angular 
correlation lengths which is reflected in the smallness of a is due to the presence of stabilizing 
extended textures within the ordering chiral field which constitutes an essential feature of 
the present investigation. 

So, if the present model with its underlying dissipative dynamics indeed provides a rea- 
sonably realistic description for the roll-down of a hot symmetric configuration into the 
spontaneously symmetry-broken T = vacuum, then we would finally conclude that the 
transient formation of DCC domains would hardly cause significant anomalous deviations 
from the Gaussian shape for the distributions of neutral-to-charged ratio fluctuations. This 



seems to be in accordance with the present experimental situation |TTJ. The model simula- 
tions, however, could lead to the expectation that the relaxation to the true vacuum after a 
chiral phase transition should manifest itself in a broadening of the width of the distribution 
as compared to standard statistical expectations, by a factor of about 1.4 — 2.4. 
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